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ABSTRACT 

Using a recently proposed approximation for multiskyrmion fields based 
on rational maps we study the masses and baryonic radii of some strange 
multibaryons within the bound state soliton model. We find the tetralambda 
binding to be stronger than previously expected. In addition, the model pre- 
dicts the existence of a "heptalambda" which is stable against strong decays. 
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In recent years there has been a continuous effort in finding minimum energy 
skyrmion solutions for increasing winding number B. A decade ago, the so-called B = 2 
torus configuration was found |]]. A few years later the lowest B = 3 and B = 4 energy 
solution were identified^]. Finally, very recently some B = 5 to 9 skyrmion configura- 
tions which are believed to represent global minima were discovered ||]. All these solutions 
could be obtained only after a sophisticated and demanding numerical work. The need 
to further investigate the physics behind these solutions has lead people to find good 
and algebraically simple approximations to them. For example, in Ref. [ffl a variational 
approximation to the B=2 torus was introduced. Even more important for our purposes, 
less than a year ago Houghton, Manton and Sutcliffe || have exploited the similarities 
between the BPS monopoles and skyrmions to propose some ansatze based on rational 
maps. They have shown that such configurations approximate very well the numerically 
found lowest energy solutions with B < 9 and conjectured they could also describe that 
with B = 17. The introduction of these approximations allows for the study of the prop- 
erties of multibaryons with strangeness. These exotics have received quite a lot attention 
since Jaffe first speculated on the possible existence of a stable H-dibaryon twenty years 
ago|§. Strange dibaryons have been extensively studied in different models including the 
Skyrme model in its various SU(3) extensions p|, |9|, [TOR . The tetralambda has been also 
considered using the SO (3) ansatz for the skyrmion field[|ll|]. In addition, experimental 
searches of these objects are in progress, particularly in heavy ion collisions at Brookhaven 
and CERN (see i.e. Refs.|L2| and references therein). The purpose of the present work is 
to extend previous investigations by considering strange multiskyrmions with B — 3 to 
9 and 17 in the bound state approach to the SU(3) Skyrme model[nj. In this approach 
strange multibaryons can be obtained by binding kaons to a multisoliton background. To 
describe such multisolitonic field we will use the rational map ansazte. 

We start with the effective action for the simple Skyrme model with an appropriate 
symmetry breaking term, expressed in terms of the SU(3)-valued chiral field U(x) as 

r = Jd A x j^Tr [dpUdrrf] + ^Tr [[U^d^U, U^d u U} 2 } j + T wz + T SB , (1) 

where F w is the pion decay constant ( = 186 MeV empirically) and e is the so-called 
Skyrme parameter. In Eq.(|]), the symmetry breaking term Tsb accounts for the different 
masses and decay constants of the pion and kaon fields while T wz is the usual Wess- 
Zumino action. Their explicit forms can be found in i.e. Ref. ||T0f . 

We proceed by introducing the Callan-Klebanov (CK) ansatz for the chiral field 



U = ^U n U K ^U„ . (2) 

In this ansatz, Uk is the field that carries the strangeness. Its form is 
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Uk = exp 
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where K is the usual kaon isodoublet K = 

^ K° 

The other component, U^, is the soliton background field. It is a direct extension 
to SU(3) of the SU(2) field u n: i.e., 

Replacing the ansatz Eq.(§) in the effective action Eq.(|l]) and expanding up to 
second order in the kaon fields we obtain the lagrangian density for the kaon-soliton 
system. In the spirit of the bound state approach this coupled system is solved by finding 
first the soliton background configuration. For this purpose we introduce the rational 
map ansatze|| 

= exp [if • 7T n F] , (5) 

with 

tt„ = | (2 $l(R n ) % + 2 %(R n ) j + (1 - |i? n | 2 ) k) , (6) 

where we have assumed that F = F(r) and R n = R n (z) is the rational map corresponding 
to winding number B = n. Here, r is the usual spherical radial coordinate whereas the 
complex variable z is related to the other two spherical coordinates (8, <fi) via stereographic 
projection, namely z = tan(6 l /2) exp(i0). 

In order to find the lowest soliton-kaon bound state we write the kaon field as|| [TXJ 



K(r,t) = k(r,t) f-7r n (z) X , (7) 

where x 1S a two-component spinor. 

Using Eqs.(|HJ^) the explicit form of the kaon-soliton effective lagrangian is 

L = J drr 2 [f n {r) tfk - h n {r) k v k' + i\ n {r){tfk - tfk) - tfk{m K 2 + V* ff )} , (8) 

where 

and 
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The term linear in time derivatives, whose coefficient is 
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K(r) = 1 + -^—^- . (10) 
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is due to the Wess-Zumino action, and 
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Here, the angular integrals I n are defined as 
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The diagonalization of the hamiltonian obtained from the effective lagrangian Eq. ( 
leads to the kaon eigenvalue equation 



-d r (r 2 h n d t 
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In our numerical calculations we use two sets of parameters. Set A corresponds to 
the case of massless pions and Set B to the case where the pion mass takes its empirical 
value m n = 138 MeV. In both cases, F n and e are adjusted so as to reproduce the 
empirical neutron and A masses. The ratio Fk/F^ and txik are taken at their empirical 
values, Fk/F^ = 1.22 and txik = 495 MeV. Once the parameters are fixed we obtain the 
soliton profiles using for each baryon number B the rational map given in Ref.[f|. Then, 
we proceed to solve the corresponding kaon eigenvalue equation. Results are shown in 
Table 1. There we list the soliton mass (per baryon number) and kaon eigenenergy for 
each baryon number for both sets of parameters. 

To get some insight on the quality of the ansatze we first concentrate on B = 2. 
In this case, the exact numerical solution []J corresponds to an axially symmetric soliton 
where the chiral angle F depends not only on r but also on 9, the angle with respect 
to the symmetry axis. A basic assumption of the rational map ansatz is, however, that 
F = F(r) for any winding number B. In addition, within that particular class of solutions, 
this ansatz proposes a very particular dependence of function of 9 and cj). For 

B = 2, the best possible solution in which the chiral angle depends only on r was found 
in Ref. j|] and applied to strange dibaryons in Ref.||10||. Comparing with the results of 
those references, we observe that M so i as predicted by the rational map ansatz is around 
20 MeV larger (for both sets of parameters) than the one in Ref. [Q Q The effect on 
the kaon eigenenergy is even smaller: the values of 62 in Table 1 are only 6 MeV larger 
than those reported in Ref.|T(J. Therefore, possible improvements upon the rational map 
ansatze are not expected to change the values of the kaon eigenenergies in a significant 
way. 

We turn now to the predictions for higher baryon numbers. As a general trend we see 
that the kaon binding energies = — e n decrease with increasing baryon number. 
However, as in the case of energy required to liberate a single B = 1 skyrmion from 

-'^As compared to the exact numerically found minima, results in Table 1 are around 50 MeV larger. 



3 



the multisoliton background p|, [§, we observe some deviation from a smooth behaviour, 
namely, Df > and Df > Dq. Consequently, such deviations will be also present in 
the multiskyrmion mass per baryon. Interestingly, this kind of phenomena has been also 
observed in some MIT bag model calculations fl4]]. There they are due to shell effects. 

Another magnitude of interest is the RMS baryonic radius of the strange multi- 
skyrmions with strangeness S. It is related to the isoscalar electric radius and gives an 
idea of the size of these exotics. It is given by 

(Vf) 2 = rL + S4, (15) 

where r 2 sol and r 2 K are the solitonic and kaonic radii, respectively, defined by 

2 r 

r 2 , = dr r 2 F' sin 2 F , (16) 

7T J 

r\ = 2 J drk 2 r 4 (f n e n + \ n ) . (17) 

The corresponding results are given in Table 2 together with the radii of the S = —1 
multiskyrmions. As expected, both the solitonic and kaonic radii increase monotonically 
as a function of B. For the radii of the S — — 1 multiskyrmions such increase is much 
slower as a result of a compensation between both contributions. Actually, for Set A the 
radius for B = 6 is even smaller than that of B = 5. This effect, however, is very likely 
to be an artifact of this parameter set since it does not appear for Set B. 

To obtain states with definite spin and isospin in the bound state model one should 
perform the SU (2) collective quantization of the soliton-meson bound system. Such quan- 



tization has to be done taking into account all the symmetries of the bound system |L0 
For B > 2, however, this is an intrincated problem which has not been fully understood 
even in the case of non-strange multiskyrmions. In this situation, we will only concentrate 
on estimating the centroid of the low-lying states assuming that, for those states, the ro- 
tational corrections are small. This approximation would become better as B increases 
due to the associated increase of the moments of inertia. Using this assumption we can 
address the issue of the stability of the strange multibaryons against strong decays. For 
the if-dibaryon we obtain 

M H - 2M A = 12 MeV . (18) 
Comparing with the results of Ref . [ 1(J we note that improvements on the description of 



the soliton background together with the inclusion of rotational effects result in an extra 
binding of ~ 40 MeV. In any case, given the various uncertainties coming from i.e. 
Casimir effects, etc., it is clear that the present calculation basically leads to the same 
qualitative result: within soliton models the H is predicted to be very close to threshold. 
The situation concerning the stability of the tetralambda with respect to the decay into 
two dilambda seems to be somewhat different, however. From Table 1 we have 

M 4A - 2M 2A = -176 MeV , (19) 
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for both sets of parameters. Therefore, it is rather strongly bound. It should be noticed 
that in this case improvements on the ansatz will decrease the binding. For example, the 
solitonic contribution in the present approximation is —200 MeV to be compared with 
the exact value —150 MeV |§. Still, the present results indicate that, within the Skyrme 
model, the binding of the tetralambda might be stronger than previously expected fTT| . 



There are some other states with vanishing hypercharge which might be stable 
against strong decays. To see that it is convenient to define the ionization energy Is as 

I B = M 1 + Mb -i - M B , (20) 

where M n is the mass of the state with B = —S = n. Within the present approximations 
M n can be calculated from Table 1 using M n = n(M so i + e n ). Results for B = 5 to 9 are 
shown in Table 3. We readily see that B = —S = 5, 8, 9 are unstable, B = —S = 6 barely 
stable while B = —S = 7 is clearly stable against one hyperon emission. It is not hard to 
check that the stability of the B = —S = 7 strange multiskyrmion is not violated by any 
other strong decay. For example, 

M 7A - (M 3A + M 4A ) = -177 MeV . (21) 

Therefore, as it stands the present model predicts the existence of a stable "heptalambda" . 

In conclusion, we have studied strange multibaryon configurations in the context of 
the bound state approach to the SU(3) Skyrme model using for the soliton background 
fields some ansatze based on rational maps. We have obtained the solutions of the kaon 
eigenvalue equation for B < 9 and B = 17. With these solutions we have calculated 
the associated baryonic radii. As expected, for S = — 1 they increase with increasing B. 
We have calculated the strange multiskyrmion masses within the adiabatic approximation. 
We have found that among the exotics with vanishing hypercharge and low baryon number 
only those with B = 4 and 7 are clearly stable against strong decays. This might be 
of interest for current experimental searches of strange exotics. Finally, we note that, 
in principle, it would be possible to study hypernuclei within the present framework. 
However, as wellknown from the case of the non-strange B = 2 system, one should 
definitely perform a quantum treatment of the 577(2) multiskyrmion before a reasonable 
description can be made. Investigations along this line are currently underway. 
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Table 1: Soliton mass per baryon unit M so i and kaon eigenenergy e (both in MeV) as 
a function of the baryon number B. Set A corresponds to F n = 129 MeV, e = 5.45, 
m n = while Set B to F w = 108 MeV, e = 4.84, m n = 138 MeV. The other parameters 
take their empirical values as it is explained in the text. 





Set A 


Set B 


B 


M vn , 

- 1 tJ - SOL 


e 


M, n , 
- 1 SOL 


e 


1 


863 


222 


864 


210 


2 


847 


244 


848 


232 


3 


830 


255 


832 


241 


4 


797 


250 


798 


238 


5 


804 


263 


808 


248 


6 


797 


267 


802 


251 


7 


776 


262 


780 


246 


8 


784 


271 


790 


252 


9 


787 


276 


796 


256 


17 


766 


284 


782 


256 



Table 2: Baryonic radii r B of the S — — 1 multibaryons as a function of the baryon 
number. The solitonic r so i and kaonic tk contributions to those radii are also displayed. 
Parameter sets A and B are defined in the caption of Table 1. All the radii are in fm. 





Set A 


Set B 


B 


r so i 


r K 


r B 


r so i 


r K 


r B 


1 


0.59 


0.44 


0.39 


0.68 


0.51 


0.44 


2 


0.81 


0.65 


0.49 


0.94 


0.75 


0.57 


3 


0.97 


0.80 


0.54 


1.12 


0.93 


0.63 


4 


1.06 


0.89 


0.57 


1.23 


1.03 


0.67 


5 


1.20 


1.03 


0.64 


1.37 


1.18 


0.68 


6 


1.29 


1.13 


0.62 


1.48 


1.29 


0.73 


7 


1.35 


1.18 


0.65 


1.55 


1.35 


0.76 


8 


1.46 


1.29 


0.68 


1.66 


1.47 


0.78 


9 


1.55 


1.39 


0.69 


1.76 


1.57 


0.80 


17 


2.05 


1.89 


0.79 


2.27 


2.09 


0.88 



7 



Table 3: Ionization energies Ib as defined in Eq.(pO|) (in MeV) for B = —S = 5 up to 9 
multibaryons. Parameter sets A and B are as in Table 1. 



B 


Set A 


Set B 


5 


-62 


-62 


6 


+36 


+36 


7 


+203 


+210 


8 


-59 


-80 


9 


-42 


-58 
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